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We investigate general braneworld models, with a non-minimally coupled phantom bulk field and 
arbitrary brane and bulk matter contents. We show that the effective dark energy of the brane- 
universe acquires a dynamical nature, as a result of the non-minimal coupling which provides a 
mechanism for an indirect "bulk-brane interaction" through gravity. For late-time cosmological 
evolution and without resorting to special ansatzes or to specific areas of the parameter space, we 
show that the —1-crossing of its equation-of-state parameter is general and can be easily achieved. 
As an example we provide a simple, but sufficiently general, approximate analytical solution, that 
presents the crossing behavior. 
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I. INTRODUCTION 

According to cosmological observations our universe is undergoing an accelerating expansion, and the transition 
to the accelerated phase has been realized in the recent cosmological past In order to explain this remarkable 
behavior, and despite the intuition that this can be achieved only through a fundamental theory of nature, we can still 
propose some paradigms for its description. Thus, we can either consider theories of modified gravity @, or introduce 
the concept of dark energy which provides the acceleration mechanism. The dynamical nature of dark energy, at least 
in an effective level, can originate from various fields, such is a canonical scalar field (quintessence) [||, a phantom 
field, that is a scalar field with a negative sign of the kinetic term Q , or the combination of quintessence and phantom 
in a unified model named quintom [f| . The advantage of this combined model is that although in quintessence the 
dark energy equation-of-state parameter remains always greater than —1 and in phantom cosmology always smaller 
than — 1, in quintom scenario it can cross —I. 

However, there are strong arguments supporting that the —1-crossing in a single, minimally-coupled scalar- field 
model, is unstable under cosmological perturbations realized on trajectories of zero measure @ . This feature, together 
with additional theoretical evidences such arc quantum-correction incorporation and renormalizability, raised the 
interest for the investigation of models where the scalar fields are non-minimally coupled to gravity [7j. 

On the other hand, brane cosmology (motivated by string/M theory), according which our Universe is a brane 
embedded in a higher dimensional spacetime [H, Q, apart from being closer to a higher-dimensional fundamental 
theory of nature, it has also great phenomenological successes and a large amount of current research heads towards 
this direction [Io| . The bulk space can contain only a cosmological constant [ll[, but string theoretical arguments 
led to the insertion of bulk scalar fields, since at the level of the low-energy 5D theory it is natural to expect the 
appearance of a dilaton-like scalar field in addition to the Einstein-Hilbert action Therefore, the cosmological 
evolution on the brane-universe is a combined effect of both the brane matter content and the bulk scalar field 
[HI, [13, HE EE [13, HH • In addition, being closer to scalar-tensor gravity theories and to the theoretical evidences 
described above, one can consider the bulk scalar field to be non-minimally coupled to the 5D curvature scalar fl9j . 

In the present work we are interested in investigating 5D braneworld models with a non-minimally coupled phan- 
tom bulk field. This scenario combines the advantages of AD phantom cosmology, together with the advantages of 
non-minimally coupled 4D or brane cosmology. In particular, the basic question is weather we can acquire an effec- 
tive cosmological evolution on the brane-universe, where the dark energy equation-of-state parameter can cross the 
phantom divide in the near cosmological past, without choosing special ansatzes for the various quantities, or moving 
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in specific, small areas of the parameter space. Since it is already known that in some special cases, non-minimally 
coupling in conventional AD cosmology can lead even single phantom- field models to experience the — 1-crossing 20] , 
it is interesting to see if this behavior can be preserved in (closer to a multi-dimensional theory of nature) brane 
cosmology. Indeed, it proves that such models not only do present the aforementioned behavior, in agreement with 
observations, but they do so for a large solution sub-class and without the need of restricting to a specific area of the 
parameter space. In particular, under the low-energy (late-time) assumptions and up to first order in terms of the 
non-minimal coupling parameter, the —1-crossing can be simply acquired if the term (fxj> on the brane- universe (where 
<j) is the phantom bulk field) changes sign. 

However, we have to make a comment about the quantum behavior of the examined model. As it is known, the 
discussion about the construction of quantum field theory of phantoms is still open in the literature. For instance in 
[2l| the authors reveal the causality and stability problems and the possible spontaneous breakdown of the vacuum 
into phantoms and conventional particles in four dimensions, arising from the energy negativity (if ones desires to 
maintain the unitarity of the theory). On the other hand, there have also been serious attempts in overco ming these 
difficulties and construct a phantom theory consistent with the basic requirements of quantum field theory [22J, with 
the phantom fields arising as an effective description. Since every warped higher dimensional model is reduced to an 
effective 4D one in low energies, the aforementioned discussion concerns the 5D phantom scenario, too. The present 
analysis is just a first approach on the —1-crossing in phantom braneworlds. Definitely the subject of quantization 
of such models is open and needs further investigation. 

The plan of the work is as follows: In section [II] we formulate general braneworld models with a non-minimally 
coupled phantom bulk field and arbitrary brane and bulk matter contents. In section IIIII we provide the general 
late-time evolution on the brane-universe and we present the main result of our investigation, namely the relation 
between the brane-universe's effective equation-of-state parameter with the values of the bulk field and its derivative 
on the brane. In section lTVl we construct a simple but quite general approximate analytical solution, which experiences 
the —1-crossing. Finally, section [V] is devoted to summarize our results. 



II. THE MODEL 



We consider a general class of single-brane models, in the presence of a non-minimally coupled phantom bulk field, 
characterized by the action: 

S= J d i xdy^Gl [ f(cb)n-A+^(Wcf>) 2 -V( ( f>)-C^+ J d±x^ {~a + c[ m) } . (1) 

Gmn describes the 5D metric while g denotes the associated induced metric on the brane, located at y = without 
loss of generality, and as usual A is the 5D cosmological constant and a is the brane-tension. The term £^ stands 

for possible forms of bulk matter apart from the phantom field, while accounts for the matter content of the 

brane-universe. For the moment we consider the coupling f(4>), to the 5D Ricci scalar 1Z, to be an arbitrary, general 
function. 

The Einstein's equations arise by variation of the action: 

f(tf>) (tZmn - IgmnTIJ - VmVjv/W + Ga/jvV 2 /(0) = ^T MNl (2) 
where Tmn is the total energy- momentum tensor: 

Tun = T^ N + + T M ^ N - G M n^ - Gn V 5 M 8 u N aS(y). (3) 
In the expression above, Tj$ N accounts for the phantom bulk field 

t mn = -VM(t>VN<f> + G MN Q (V0) 2 - V(<t>)) . (4) 

Similarly, T^ N stands for the bulk part of the energy-momentum tensor, while accounts for the brane-matter 

content. Finally, varying the action in terms of the phantom field, we obtain its evolution equation: 

o dV _ df da „ , , 

- v - -JT + n yi - -ji s (y) = °- ( 5 ) 
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For the metric we will use the following form [2c 



ds 2 = -n 2 (y,t)dt 2 + a 2 {y,t)~/i 3 dx' l dx j + b 2 (y,t)dy 2 



(6) 



which corresponds to a maximally symmetric, induced Friedmann-Robertson- Walker geometry on the brane. Thus, 
the isomctry assumption along three dimensional x-slices, including the brane, allows as to consider the phantom bulk 
field depending only on the fifth coordinate. Furthermore, using square brackets to denote the jump of any quantity 
across the brane ([Q] = Q(y+) — Q(y_)), and assuming S 1 /Z2 symmetry across it, we restrict our interest only in the 
[0, +00) interval, obtaining 



[Q 1 ] = 2Q'(0+). 



(7) 



This relation is going to be used in the elaboration of the discontinuities of the derivatives of various quantities at the 
location of the brane [24[ . 

In these coordinates, the phantom energy-momentum tensor writes: 



(<t>) 



MN 



%?r+n 2 v 













^</> 2 + - v 



2b 2 ' 



(8) 



b" 

TP 



b 2 V 



where primes and dots denote derivatives with respect to y and t respectively. For the brane-universe we assume 
that it contains a perfect fluid with equation of state p = w m p, where p and p are its energy density and pressure 
respectively. Thus, the brane energy-momentum tensor reads: 



T 



MN 



pn 2 
pa 2 ^ 




(9) 



Similarly, assuming an ideal fluid for the arbitrary forms of bulk matter, with energy density ps and pressures Pb 
and Pb (since the pressure on the fifth dimension can be different), we acquire: 



-MN 



pBn —n 2 P 5 

P B a 2 ll0 _0 
-n 2 P 5 P B b 2 



(10) 



Note that we have allowed for an energy-exchange function P 5 between the bulk and the brane-universe [16[ . Finally, 
all quantities are considered as functions of y and t, apart from the brane-fluid's p and p which depend only on time. 

In order to focus on the cosmological evolution on the brane we use the Gaussian normal coordinates (b(y,t) = 1) 
[l4| . Thus, the non-trivial five-dimensional Einstein equations consist of three dynamical: 
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and two constraint equations 
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Note that we have allowed for a ^-dependence of the brane tension <j(<f>). Finally, as usual the 5D Ricci scalar 1Z is 
given by 



K = 3- 
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In order to obtain the boundary conditions for the aforementioned cosmological system, we integrate the 00 and 
ii components of the 5D Einstein equations around the brane, making use of ([7]). Thus, we result to the following 
junction (Israel) conditions (setting also n(0,t) = 1 without loss of generality): 



-6 A 

a 



d0 



4 A + 2n' f - 2 f 
ao d<p 



20o + 4 |r 



>>'o = ^(p-cr) 



da 



3^ ) = - — 

ao 



(17) 
(18) 
(19) 



where the index denotes the values of the corresponding quantities at the location of the brane. For notation 
simplification in the following we omit it, i.e a, n and <f> and their derivatives, stand for the corresponding quantities 



on the brane. Furthermore, we call and and f a the terms ^ 
Finally, equations (|17| can be re-written as: 



and 4? 



respectively. 



a' = l (<r + 2/ / g ) l (3p-p) 1 (p + p) 
a 2 u 8 u 16 / : 



(20) 



, l(a + 2Uf a ) l(3p-p) , 3 (p + p) 



2 u 



u 16 / 



(21) 



b' = -U 



(3p-p) (3fU-4aU) 



(22) 



where we have defined 



u = 6/ + 16/| 



(23) 



III. GENERAL LATE-TIME COSMOLOGICAL EVOLUTION ON THE BRANE-UNIVERSE 

In this section we develop the formalism for the late-time brane-evolution investigation, following [l5|, [Q. First of 
all, using the boundary conditions (Ti"T)) - (fll))) , together with (fT3j) at y = 0, we acquire the energy conservation equation 
on the brane: 



p + 3-(p + p) + 2P 5 = 0, 



(24) 



which in the case P5 ^ describes the direct bulk-brane energy flow, which can have both signs. Similarly, using 
equation (jTTJ) at y — 0, we obtain the cosmological evolution on the brane: 



3/ 
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(25) 
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Finally, the use of (fl3]) at y — provides the phantom field evolution on the brane: 

T 



dV 
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2^ 



a 



+<f>' [ n' + 3- ) +6/ 
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y-2/^(3— +fi"). (26) 



The terms p'(t),a''(t),n''(t) denote the unknown, "non-distributional" parts of the corresponding derivatives [18l |. 
Although they could be set to zero, we consider them to depend on the structure of the bulk, a" and n" can be 
expressed in terms of cj>" and the standard quantities, i.e. a(t),<p(t), their time-derivatives and the various matter 
densities. Assuming f a = and following |18|, the corresponding elimination leads to the evolution equations on the 
brane: 
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(27) 
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where 



da 



and 



d 2 f 



We mention that in order to simplify the above two equations we have focused on 



the late-time approximation, namely neglecting y0 2 -terms [la , |25| since late time is equivalent to low-energy. Finally, 
note that the phantom evolution equation still contains the unknown function <p" . In order to provide a specific 
analytical solution sub-class we have to use an ansatz for it, and this will be done in the simple example of the next 
section. 

Equations (|24|) . (|27|) and (f28|) describe the cosmological evolution of a brane- universe with arbitrary curvature, in 
the case of arbitrary bulk and brane contents, for a general coupling f(<fi), in the low-energy approximation. In order 
to proceed we have to make a choice for the form of f ((/))■ We choose the well-studied quadratic form: 



Ml\l 



(29) 



where £ is the coupling parameter [l9j , and M5 the 5D Planck mass. This ansatz is a good approximation to a general 
coupling function for small phantom field values. Using (j2"9"j) we obtain: = — M|£0 and u = 6M| ^1 — §0 2 ^ + 
16Mf£ 2 2 . 

As a next step we proceed to the dark energy field formulation [l5| . introducing the dark energy field %(t) in a way 
that the second-order equation (|2"T|) is replaced by two first-order ones: 



a \ k 
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where (3, 7 and A are given from 
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A simple verification of the above formulation can be obtained by differentiating (|30p . substituting in (|3ip and using 
In this way the initial equation (|27p can be recovered. 
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Under the constructed formulation, equation (|30|) corresponds to the conventional AD Friedmann equation of the 
(brane)-universe, and the function A is just the effective AD cosmological constant of the brane. In the minimal 

coupling case, i.e when £ = 0, we obtain A = 12 * /:i + 12 ° fl/a j = const (i.e. A = 0) and thus the universe's dark 

energy behaves like a cosmological constant. In addition, one can fine-tune A to zero [9| as: 

A =^( A+ i!y =°< ^ 

which corresponds to the familiar Randall- Sundrum fine-tuning and to the known brane solution with no dark energy. 

In the non-minimally coupling case we observe that A is not a constant anymore but it has acquired a dynamical 
nature (A ^ 0) due to the dynamical nature of the coupling /(</>(£)). Therefore, A corresponds to the effective AD dark 
energy of the brane-universe. Surprisingly enough one can acquire an analytical expression for its equation-of-state 
parameter w e ff, up to second order corrections in terms of the coupling £. In particular, differentiating (|34p , using 
the ansatz (1291) and keeping only first order £-terms we obtain: 

A - i(t>4>\ = => A + 3- ( - IcHt I A = 0. (36) 

a \ 3 a J 

This relation is independent of A and £ normalizations (since some authors consider 2M| instead of Af|) and thus 
general. Therefore, we can straightforwardly make the identification: 

Weff = -1 - g£#| ■ ( 37 ) 

Relation (|37[) is the main result of the present work, and provides the equation-of-state parameter for the dark 
energy of a brane-universe, embedded into a bulk with a non-minimally coupled phantom field and arbitrary bulk 
and brane matter contents, up to first order in terms of the coupling £. In the minimally-coupled case (£ = 0) 
we acquire w e ff = — 1, that is we verify that the dark energy is simply just a cosmological constant. However, in 
the non-minimally coupled model, which is the case of interest of the present work, we observe that the obtained 
cosmological behavior is very interesting. In particular, if (fxp is negative, then — 1 < w e /f, that is the universe's dark 
energy behaves like quintessence. It is interesting to notice that we obtain an effective AD quintessence behavior, 
although we have the presence of a phantom bulk field. On the other hand, if <jxp is positive, then w e ft < — I and 
dark energy behaves like a AD phantom. 

The most interesting case is when <jxj) changes sign during the cosmological evolution. In this scenario w e ff crosses 
the phantom divide. In addition, if initially we have <jxf> < and after some particular moment we have <jxf> > 0, then 
the — f-crossing can take place from above to below, that is consistently with observations. We mention that the 
crossing can be achieved with the use of only one bulk field, thus the constructed model is quite economical. This 
feature was already known to hold in AD single-field, non-minimally coupled cosmology, but only in specific cases j20| . 
However, the non-trivial fact that it is maintained in the higher-dimensional, and thus closer to a fundamental theory 
of nature, brane cosmology, and without the need of restricting to specific, small areas of the parameter space, provides 
an additional argument for the significance of the non-minimal coupling for the description of nature. Finally, note 
that the dynamical nature of the AD effective dark energy and the possible brane-universe acceleration, is obtained 
independently from the possible direct energy flow between the bulk and the brane. In other words, the non-minimal 
coupling provides a mechanism for an indirect "bulk-brane interaction" , through gravity. 

In conclusion, in this section we have formulated the general late-time (low-energy) cosmological evolution on the 
brane-universe. In order to acquire a specific example one can either solve the aforementioned equations numerically 
or choose some simple ansatzes for the scale factor, the bulk field and the matter terms, and solve the equations 
analytically. In the next section, we use some general but simple ansatzes, in order to present more transparently the 
cosmological implications of the constructed model. 

IV. CROSSING 1: A SIMPLE APPROXIMATE SOLUTION 

Let us now present a simple, but sufficiently general, approximate, analytical cosmological solution. Since we focus 
our analysis at late times, it is reasonable to consider the widely-used power-law ansatzes. Thus, for the bulk field we 
assume: 



(38) 
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with C Q , Op being constants. Similarly, for the scale factor we can safely consider the solution: 

a(t) = C 3 t Ul . (39) 

Now we have to choose the exchange term P 5 that is present in as well as the bulk pressure Pb that appears 
in (|25p . which arc both functions of time, corresponding to the values of Ps(y, i) and P§(y,t) on the brane. The 
energy-momentum conservation VmTjS = cannot fully determine Pb and P5 and a particular model of the bulk 
matter is required. Although we could consider the ansatz [l6[ (see also [Hj]): 



P 5 (t) = C 5 



a(t) 



a(ty 



(40) 



for the investigation of this section and for simplification of the results it is adequate to set P 5 = 0. For the bulk 
pressure we assume Pb = const, and the constant can be absorbed in a re-definition of the bulk potential. This can 
be chosen as (l7| : 



V(<P) = 



(41) 



which is a general form, consistent with the brane stabilization mechanism 27J ■ Finally, for the matter content of the 
universe we have already assumed an ideal fluid, with energy density and pressure connected by: 

P = w m p. (42) 

Thus, using (14"2")) we can easily solve (|2~4"|) obtaining: 



C„ 



P = 



(43) 



a 3(l+1« m ) ' 

with a(t) given by (|3Uj) . 

We mention that the considered ansatzes for a(t), P5, Pb and V(<j>) are not crucial for the obtained results, and 
they are chosen just for presentation reasons. The only necessary behavior is a (b(t) with a derivative that changes 
sign during the evolution, as we have discussed in the previous section. The ansatz (|38[) can fulfill this conditions 
in a simple way, but any other ansatz with a derivative sign-flip could be equivalently used. Finally, if one use an 
arbitrary ansatz for cf> without a derivative sign-flip, then according to (|37p he will acquire a w e tf lying always on the 
same side of the phantom divide during cosmological evolution. 

The final step is the substitution of the aforementioned ansatzes, together with (|43|) . into the cosmological equations 
([30|) and (f3Tj) . In order to simplify the calculations we consider only the flat brane- universe case (k = 0), and since 
we are investigating the late-time behavior we only keep terms up to 0(t~ 2a ) and 0(t~~ 213 ) [18]. Finally, as stated 
in section UTT1 we keep terms up to 0(£ 2 ). Therefore, within these approximations, the two cosmological equations 
become: 
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+ O{t- 2a ) + O{t- 2(3 ) + O(£, 2 ) = 0. (45) 



A last assumption concerns the unknown function 4>" . Although its specific form is not important, in order to simplify 
the calculations and following [l8[ we consider: 
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^ lcx+2 f-13+2 ' 

In order for equation (|4"4"|) to be satisfied for all t we require: 
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Note that (|48[) is the standard Randall- Sundrum fine-tuning for a vanishing cosmological constant on the brane, in 
the minimally-coupling case in the absence of matter. Under these conditions, the requirement of satisfaction of (|45[) 
at all times leads to: 



fr 2 
12M§ 



C a a(a + 1 - 2,vx) - 16M 5 3 C* Q v x (2v x - l)f - C 0Q = 
Cp /3(/3 + 1 - 3^) - WMlCp v x {2vx - \)i - C H3 = 0. 



(50) 

(51) 
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Conditions (147l) - ([52|) are necessary in order for the selected ansatzes to form a self-consistent cosmological solution. 
The quadratic equations ([ST]) and |52|) lead to: 
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As we can see, if v\ > 1/2, that is according to (j4"7| if w m < 1/3, the exponents a and f3 are real for any value of £. 
On the other hand for v\ < 1/2 ( i.e iu m > 1/3) the condition a,/? € R leads to a restriction in ^-values, namely: 



£ < rnin(£i,£ 2 ), 



where 



6 =Cc 
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In these expressions we have used the definition £ c = 



— 3 



16M? 



namely the conformal value of the non-minimal coupling 
> 1/3, the value of £ must be sufficiently small. 



parameter £. Therefore, in the (non-conventional) case w 

Let us now examine the cosmological behavior of the simple, approximate, but quite general, model of this section 
According to ([57)l . and using the ansatzes (13"5|) we obtain: 



"eff 



= -1 



1 

3 v\ 



aC a 
t a 



(58) 



with a and /3 given by (fS"3"|) and ([54^1 . Note that due to the parameter eliminations through conditions (|4"71) - (f5"2")) . one 
can either impose a value for w m and acquire the value of C p , or the opposite. Since the direct observable is w m we 
prefer to use its value as an input, and thus eliminate C p from Wetf's calculation. 

As we have already mentioned, for £ = dark energy behaves like a cosmological constant (w e ff — — 1 = const), 
but the non-minimal coupling gives it a dynamical nature. Assuming <j){i) > (i.e. choosing C ai Cp > 0), then if 
4>{t) changes sign during the cosmological evolution then w e ff crosses the phantom divide —1, while if <p(t) preserves 
the same sign, w e tf lies always on one side of the phantom divide. However, one could consider cases where <p(t) 
changes sign or where both <fi(t) and <j>(t) do. Finally, as can be seen from (|58[) and (|53|) . (|54p the decisive role for the 
determination of w e ff is played by the ^-parameters (C a ,C^), w m and of course £, with the rest parameters being 
non-relevant (for realistic choices). 

In order to acquire a better comparison with the observational data, in the following we provide the evolution of 
Weff versus the redshift z, given by (1 + z) = ao/a(t) with ao the present value, where we fix C3 in order to acquire 
do = 1. Furthermore, for the equation-of-state parameter of the matter content of the universe we use the value 
w m — 1/3, since a relativistic matter is the most natural choice. This value for w m has an additional advantage, 
namely it leads to 1/1 = 1/2 (see (|47|) ) and thus A/5 is eliminated from the final results (see (|53|) and (|54[) ) . Therefore, 
since M5 defines the units of the present model, its disappearance from the final relations ( (|4T)l . (|5"3"|) , (f5"4")l ) and ([55]) ) 
allows us to choose the units at will. Thus, we choose the units in order for the parameters C a , C^a, Cp, C^p to be of 
order 1. 
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FIG. 1: (Color online) w e ff versus z, for w m 
0, Cj,p = 0, while the lower one to C a = 0, C^a 



■■ 1/3 and £ = 0.1. The upper curve corresponds to C a 

0, Cj3 = 1, = 1. 
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In fig. Q] we depict 



U eff 



versus z for £ = 0.1, and two choices of the parameter four-plet C Q , C^ Q , Ca, C^a. This 



value for £ is consistent with 0(£ 2 )-calculations, while the CVparameter choices correspond to keeping only one term 
in <j){t), thus its derivative preserves the same sign and therefore w e ff behaves like in AD quintessence (upper curve) 
or in 4D phantom (lower curve) paradigms. 

In fig. [2] we depict w e ff versus z for four choices of the parameter-group £, C a , C^ a , Cp, C$p. In this case we are 
interested in acquiring a sign-flip of <j), and thus without loss of generality we keep the plus sign in a-solution in (|53|) 
and the minus sign in /^-solution in (|54p . As we observe, w e ff crosses the phantom divide from above to below in the 
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FIG. 2: (Color online) w e ff versus z, for w m = 1/3. The values of the parameters £, (7 a , C^q, (7/3, (7^,3 are shown in the inset. 



recent cosmological past, as required by observations. This behavior is qualitatively independent of the values of the 
parameters £, C a , C^a, Cp, C^p, however the precise values of z c (redshift at the crossing) and of w e //o (present value 
of w e ff) do depend on them. Thus, we see that in the simple, but quite general, example of the present section, the 
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— 1-crossing can be achieved relatively easily, without a restriction to a small area of the parameter space. This feature 
makes braneworld models with a non-minimally coupled phantom bulk field a good candidate for the description of 
the current universe acceleration. 



V. CONCLUSIONS 



In this work we examine general braneworld models, with a non-minimally coupled phantom bulk field and arbitrary 
brane and bulk matter contents. Imposing the low-energy, that is late-time assumptions, and performing the calcula- 
tions up to first order in the non-minimal coupling £, we provide a general relation that connects the equation-of-state 
parameter w e ff of the AD effective dark energy, with the values of the phantom field <p and its derivative at the 
location of the brane-universe. For £ = 0, w e ff is always —1 and the dark energy of the brane-universe behaves like 
a cosmological constant, as expected. However, when the non-minimal coupling is switched on, the brane-universe's 
dark energy acquires a dynamical nature. 

In particular w e ff is related to £, to the scale factor and its derivative, and to 4> and its derivative on the brane 
(relation (|37[) h Thus, if 4>4> preserves the same sign during the cosmological evolution, then w e ff remains always 
on the same side of the phantom divide and the universe's dark energy behaves like a quintessence or conventional 
phantom. On the other hand, if <p4> experiences a sign-flip at some particular time, then w e // crosses —1, in agreement 
with observations. This behavior is general and, up to first order in £, it is independent of the bulk and brane matter 
contents. That is, the non-minimal coupling provides a mechanism for an indirect "bulk-brane interaction" , through 
gravity. Furthermore, the crossing behavior appears without the need of special ansatzes for the various quantities, 
or of the restriction to specific, small areas of the parameter space. 

In the present model, the — 1-crossing arises from a single, non-minimally coupled phantom bulk field, and thus 
it is more economical than the multi-field models of conventional cosmology. This behavior was already known to 
hold in AD single-field, non-minimally coupled cosmology, but only for specific cases [20j. The fact that it is not only 
maintained in the higher-dimensional (and thus closer to a fundamental theory of nature) brane cosmology, but it 
appears for a much larger solution sub-class and parameter sub-space, is definitely a novel feature and an advantage 
of the model. This result provides an additional argument for the significance of the non-minimal coupling between 
gravity and scalar fields, for the description of nature. However, the possible quantization difficulties of such phantom 
scenarios is an open problem and needs further investigation. 

Acknowledgements: One of us (E.N.S) wishes to thank Institut de Physique Theorique, CEA, for the hospitality 
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